AN ALGORITHM FOR COMPUTING THE CENTERED 
HAUSDORFF MEASURES OF SELF-SIMILAR SETS 

MARTA LLORENTE AND MANUEL MORAN 

f"*) Abstract. We provide an algorithm for computing the centered 

P\j HausdorfF measures of self-similar sets satisfying the strong sepa- 

1— H ration condition. We prove the convergence of the algorithm and 

^ test its utility on some examples. 
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1. Introduction 



We present an algorithm that takes as input a hst \I/ = {fi, f2, ..., fm} 
of contracting simihtudes in M" satisfying the strong separation condi- 
tion (see section[2]) and gives as output an estimate of the s-dimensional 
centered Hausdorff measure, C^{E), of the self-similar set E generated 
by ^. Here s is both the similarity and Hausdorff dimension of E 



CO and it can be computed from the contracting factors of the similitudes 



t^^ (see ( 2.1[ )). To our knowledge this is the first attempt at automatic 

O computation of the exact value of a metric measure (e.g. Hausdorff, 

,__! packing, spherical, centered, ...), a topic which has generated a consid- 

^ erable quantity of research ([1], [1], [7], |8], [12], [H], [H], [19] [20], [21], 

[22], [23], etc.) 
The centered Hausdorff measure is a variant of the Hausdorff mea- 
C^ sure. These measures differ mainly in the natures of the coverings used 

in their definitions. In the case of the centered Hausdorff measure we 
consider only covers by closed balls B{xi) centered at points Xi in the 
given set. 
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The standard definition of C^ on subsets of M" consists of two steps. 
Given A C M", we first compute the premeasure Cq as 
(1.1) 

CqA = hminf \ ^ \Bi{xi)f : A C UBi and \Bi\ <6,Xi e A for all i 

Because the suppression of good candidates for the Xj's may cause an 
increase of the infimum, this premeasure is not monotone, although it 
is cr— additive (see [15]). In order to avoid this difficulty we define 

(1.2) C'A = sup {C'^B : B CA,B closed} . 

The set function C^ so obtained is a metric measure. It turns out that 
the centered Hausdorff measure is bounded by constant multiples of 
the ordinary Hausdorff measure (see [E]). More precisely 

2-'C'{E) < H'{E) < C'{E) 

and so the centered Hausdorff dimension and the ordinary Hausdorff 
dimension coincide. In particular, for the self-similar set E, we have 
that < C^{E) < oo. Thus C"* is a nice measure, but there remains 
the question of how to compute C^A for some subset A G E. A 
simplification comes from the following observation: for any Borel set 
A C M", 

C'A = C'EfiA, 

where is /x the so called natural or empirical probability measure on E. 
Therefore, the problem of computing C^ on E reduces to computing 
C^E. Observe that the yU— measure of any open subset A G E with 
boundary dA having null /i— measure (and in particular of any open 
ball) can be easily obtained with arbitrary accuracy and, hence, the 
/i-measures of compact subsets of E with /x— null boundaries can also 
be computed through their complements. This gives a vast class C of 
Borel sets with computable /z— measure. In turn, the C*— measure of 
any C*— measurable set can be approximated with arbitrary accuracy 
by the C'^— measures of closed sets (see [6] Theorem 1.6 b). 

How, then, to compute C^El Given the definitions, the task seems 
out of reach. The first obstruction comes from the need of the second 



step 1.2 in the definition of C*. However, in [TU] it is proved that, 
for any subset A of a self-similar set E as above the measure and the 
premeasure coincide: 

1.3. Theorem (Theorem 3 in fTO]). Let A be either a closed or an open 
subset of a self-similar set E satisfying the open set condition. Then 
CqM = C'A. 
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With this resuh available, C* seems easier to compute than the Haus- 
dorff or spherical Hausdorff measure. Namely, the differences between 
these three measures are that, for the Hausdorff measure, one optimizes 
among coverings by arbitrary convex sets; for the spherical Hausdorff 
measure one uses arbitrary balls; and for C^ one uses coverings by balls 
with centers in E. For the Hausdorff and the spherical Hausdorff mea- 
sures, the classes of available coverings are larger, and hence it is more 
difficult to find optimal coverings. 

The second step which permits the computation of C^ was also given 
in [To]. There it is proved that computing C* is equivalent to finding 
a centered ball with optimal inverse density: 

1.4. Theorem (Theorem 5 in [lOj). Suppose the invariant set E of 
the system ^ satisfies the open set condition, with dim// E = s and 
\E\ = R, and let fi be the normalized Hausdorff measure on E. Then 

Moreover, if \1' satisfies the SSC then 

(1.6) C'E = min I } , „, : x e E and c < d < R 

{n{B{x,d)) 

where c := minjjgM,i^i(/i-E, fjE) and R := \E\. 

From now on, B{x,d) denotes the closed radius d ball centered at 
X eM". 

The statement (|1.6) is crucial in our approach since it gives that if 



(2dr 



B{x,d) is a ball of maximal inverse density, then C'^{E) — (Btxd))- 
Therefore, to obtain the value of C^{E) we need only to find an op- 
timal ball and compute its density. This is precisely how the algo- 
rithm proceeds. It searches for balls that maximize the density func- 
tion hs{x,d) = ,L j ^.. . Actually, by the self-similar tiling principle 
(see [13]), we know that finding an optimal ball is equivalent to finding 
an optimal covering. This is so because we can get an optimal covering 
tiling the set E with balls of optimal density. 

E may be tiled, without loss of /^-measure, by tiles similar to a given 
tile B. By similar we mean that the tile is an image of B under a 
composition of similitudes in \E'. The only condition to be imposed on 
B is that it be closed and have /i-B > 0. 

Our algorithm computes C^E and provides an (approximate) ball 
of maximal density, together with an optimal covering of E by balls 
centered at E. 
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Now we describe the main steps of the algorithm. This is done 
rigorously in section [3j Recall that the goal is to find the ma ximal 



value of hs{x,d) = ,L j ^.. for x G -E and d G [c, R] (see (1.6)). At 
step k, the set E is replaced by a finite set of points {A^} such that 
UAk = E and the measure /i is replaced by a discrete probability 



measure Hk supported on Ak and converging weakly to /x (see (3.4) 



(3.9) and Lemma 4.1 (iv)). The objective now is to find the maximum 



of the discrete density function 

7/7/ \\ {'^d{x,y)Y . , . 

hk{x, d{x, y)) = — -— — r^ with x,y e Ak. 

fik{B{x,d{x,y)) 

Here d{ , ) stands for the Euchdean distance. For each x E Ak the algo- 
rithm searches for the maximal value of hk{x, d{x, y)) for y G Ak- Once 
this has been found for every x E Ak, the algorithm finds the maximum 
of these values with respect to x. Thus we only need to compute ex- 
actly fik{B{x, d{x, y)) for every y G Ak- To this end, the points y E Ak 
are listed in order of increasing distance to x, and thus the points pre- 
ceding a given point y in the list always belong to the ball B{x, d{x, y)). 
It is not hard to see that the exact value of fik{B{x, d{x, y)) is obtained 



from the place of y in this list (see (3.5) for the homogeneous case 



and (3.11) for the general case). It remains to show that in the limit 
k /^ oo this process converges to C^{E). This is done in section |4J 
The convergence is shown in two steps. First, by means of the Markov 
operator M associated to the set E (see section [2]for notation and defi- 
nitions), the measures Hk are shown to converge weakly to the invariant 
measure /i. The basic properties of these measures yield a sequence of 
pairs of points {xk, yk) G A^ x A^ such that hk{xk, d{xk, yk)) — ^ C'^{E). 
However, there is no reason why these {xk,yk) should optimize hk. 
Nonetheless, we are able to show that this holds asymptotically, which 
is enough for our purposes. An interesting technical point is that in the 
proofs an essential role is played by a result of Mattila [TT] implying 
that fi{dB{x,d)) = 0. 

In section |5] we apply the algorithm to treat several sets whose cen- 
tered Hausdorff measures were available in the literature. It is remark- 
able that in all these cases the optimal value (and also optimal ball and 
covering) is attained at an early iteration. The algorithm also yields 
conjectural values (which, in many cases, are upper bounds) for sets 
whose measure is unknown. Research is in progress to explore the rate 
of convergence and show that the method yields the precise values of 
C'^{E). Preliminary results seem to indicate that, for self-similar sets 
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with less than four similarities, four decimal digits of accuracy are at- 
tainable by personal computers without any serious effort to optimize 
the code's design. 

2. Preliminaries 



Let \1' = {/i, /2, •••, /m} be a hst of contracting similitudes on M", 
\fi{x) — fi{y)\ = ri\x — y\ where < rj < 1. The unique non-empty 
compact set satisfying 

where S'\E'(X) = |J fi{X), is called the self-similar set associated to 

i<=M 

^. Sometimes we shall refer to E as the attractor or invariant set of 
the iterated function system (IFS) \E' = {/i, /2, ..., /„}. 

We shall use the following notation. Let M := {1, ...,m} and 

M^ = {ifc = (^1, ..., tk) -.tjeM Vj = 1, ...m}. 

For ifc = iii2--ik G M'^, we write 

/ife = fh° fi2° ■■■° fik^ 

ty , rtri ty ly 

' Ife — ' «l' «2---' *fc' 

and for A C M", we write 

A, = fiM)- 

We shall refer to the sets Ei^^ = f\,.{E) as the cylinder sets of generation 
k. 

Throughout the paper we shall assume that the system \E' satisfies 
the strong separation condition (SSC). That is, fi{E)nfj{E) = for all 
i 7^ j, where i,j G {1, ...,m}. The Hausdorff dimension of E, dim^^ii^, 
is given by the unique real number s such that 

m 

(2.1) E^^^ = i- 

Moreover, the Hausdorff measure of E is finite and positive. 

We shall denote by fi the natural probability measure, or normalized 
Hausdorff' measure, defined on the ring of cylinder sets by 

(2.2) ^(Ei) = r? , 

and then extended to Borel subsets of E. The measure n is regular in 
the sense that there are positive numbers a and b such that 

(2.3) ar'' < n{B{x, r)) < br'' for x E E and < r < 1. 

Let P(]R") be the space of probability measures on M". The well 
known fact that the measure /x is the unique invariant measure for the 
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Markov operator (see, for example, [2]), plays an important role in our 
proofs. Let M : P(]R"') — )■ P(]R") be the Markov operator associated 
with the IFS \1/ with probabilities {r^, ...,r^}, 

(2.4) M(z/) = rfz/o/-i + ...+r>o/„i Vz/ G P(M"). 

Then fj, is the unique probability measure satisfying 

3. Description of the algorithm 

In this section we introduce an algorithm to compute the centered 
Hausdorff measure for self-similar sets satisfying the SSC 

Given A C M", \A\ stands for the diameter of A. 

In [TU] it is proved that if ii^ is a self-similar set satisfying the SSC, 
then 

(3.1) C'E = min I ,i, ' „, : x e E and c < d < R 

[fi{B{x,d)) 



where c := inmij^M,ij^j dist{fiE, fjE) and R := \E\ (see Theorem 1.4). 



Our method depends strongly on (3.1) as, to find the value of C^E, 
we construct an algorithm for minimizing the value of 

h ( r^\ - (2^)' 

'^'^''' ' ■' fx{B{x,d)) 

when X & E and c < d < R. 

The idea is to construct a sequence {Ak} of countable sets and a 
sequence {fik} of discrete measures supported on the A^ such that 
the /ifc converge weakly to fi and U^j^A^ = E, where A stands for 
the closure of A. This allows us to construct another sequence {rhk} 
converging to C^E by choosing on the kth step a pair {xk, yk) & AkxAk 
satisfying 

rfifc := hk{xk,dist{xk,yk)) = min{hk{xk,dist{xk,yk)) ■ {xk,yk) e AfcxAfc}, 

where /ifc(x,rf):=^^^^|||^. 

We describe first the algorithm for self-similar sets where all the 
contraction ratios coincide, as this case illustrates better the central 
idea of the construction. After this we shall explain the modifications 
needed to treat the case of unequal contraction ratios. Observe that if 
Ti = Tj := r for all i ^ j, the invariant measure /i satisfies that 

(3.2) ^(E,J = r''^ = —^ Vi.GM^ 

3.3. Algorithm. (Homogeneous case: Tj = Vj := r V? 7^ j, i,j = 1, ...,m) 
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1. Construction of A^. Let Ai = {xi E M" : fi{xi) = Xi, 
i = 1, ..., m} be the set of the fixed points of the similitudes in 
\1/. For k e N+, let At = S'^^Ak-i) be the set of m^ points 
obtained by applying S'\E'(a;) = |J fi{x) to each of the m''"^ 

points of Ak-i- 

2. Construction of /i^. For all k E N+, set 



(3.4) /ifc(x) = — \/xeA 



m 



k. 



Thus, 



j=i 
is a probability measure with spt{fik) = A^ = {xi, ...,x„fc}. 
3. Construction of nik 

3.1 Given xi G Ak, compute the m^ distances dist{xi,x) for 
every x G Ak.. 

3.2 Arrange the distances in increasing order. 

3.3 Let di < d2 < .... < d^^k be the list of ordered distances. 
Since for each j = 1, ..., m'', B{xi, dj) contains j + t points 
of Ak, where t is the cardinality oi Ak r\ d B{xi,dj), i.e. 
di = di+i = ... = di+t 7^ di+t+i, there holds 

(3.5) fik{B{xi,dj)) = ^^. 

3.4 Let ijfc(/) = {ii{l), ...,ik{l)) G M^ be the unique sequence of 
length k such that Xi = fif.(i){x) for some x E Ai. We shall 
use the notation a;ij.(i) := f-i^(j_){x). Compute 



(3.6) 



{2djy _ {2djy _ m'' {2dj 
fik{B{xi,d,))~^r~ l + t 



.k-li 



only for those m {m — 1) distances dj satisfying the fol- 
lowing condition. 
Condition: If 

dj = dist{xi,y) 

for some y E Ak such that y = fjf.{z) = fj^...j^{z) for some 
z E Ai, then 

(3.7) n(O^Ji- 

3.5 Find the minimum of the m^~^{m — 1) values computed in 
step 3.4. 

3.6 Repeat steps 3.1-3.5 for each xi E Ak, I = 1, ...,m^. 
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3.7 Take the minimum of the m^ values computed in step 3.6. 

3.8. Algorithm (General case). The main difference between this case 
and the previous one is that the values of the measures fik are different. 
The structure of the algorithm is the same in both cases. Consequently, 
we shall list only the changes needed to find the measure when the 
contraction ratios are unequal. 



1. In step 2, replace (3.4) with 

(3.9) f^k{x)=rl \fxeAk, 

where i^ = (ii, ...,ik) G M^ is the unique sequence of length k 
such that X = fif,{y) for some y E Ai. Thus, /i^ is a probability 
measure with spt{fik) = ^k = {xi, ..-^Xmk}. If, for every j = 
1, ..., m*^, we denote by ifc(j) G M'^ the unique sequence of length 
k such that Xj = fi^{j){y) for some y G Ai, then we can write 



(3.10) 



fJ'k 



Z^^ik{j)^^r 



i=i 



2. In step 3.3, replace (3.5) with 



j+t 



(3.11) 



IJ,kiB{xi,dj)) = ^■ 



ifc(«<j)' 



9=1 



3. 

(3.12) 



gj G M" is such that dg 



where iki^ 

t. 

In step 3.4 , replace (3.6) with 

{2d,y 



d{xi,Xi)ioTal\q = l,...,j + 



{2d, 



lik{B{xudj)) Y, 



9=1 h{iq) 



3.13. Remark. Note that, in (3.6) and (3.12), we only compute the 
values of the inverse density function for those distances between points 
that belong to different basic cylinder sets. Namely, if dj = dist{xi, Xp) 

for some Xp G A^, then '^ /2dY i^ ^^^y computed by the algorithm 
when xi G Ei and Xp G Ej with i ^ j. 

3.14. Notation. In the rest of the paper we shall use the following 
notation. Let A^ = S'\E'(Afc_i) be the set of m'' points obtained after k 
iterations with Ai = {x G M" : fi{x) = x,i = 1, ...,m}. We write 



A:-- 



^T=iAk. 
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For k G N^ and for each xi G Ak, let D[ be the set of m^'^im — 1] 



distances satisfying condition (3.7) in the construction of the algorithm, 
denote by 

the set of the im?^~'^{m — 1) distances that appear in step 3.5, and write 

D := U- oA. 



Observe that D only takes values in the interval [c,R] (see (3.1)). 



From now on, we shall assume, without lost of generality, that R : = 
\E\ = 1. 

4. Convergence of the algorithm 

4.1. Preliminary results. The next two lemmas collect some basic 
results needed in the proofs of our theorems. We shall prove only those 
statements that do not follow directly from the construction of the 
algorithm. 

4.1. Lemma. 

(i) For every x & E there exists a sequence {xk} C Ak such that 

Ymvk^^Xk = X. 
(ii) For every k G N+, 

Ak C Ak+i- 

(iii) Let k G N+, x G Ak, and i^ G M^ be such that fi,^{y) = x for 
some y & Ai. Then 

(4.2) f,k{x)=fik{Ei,) = fi{E-J. 
(iv) For every k G N^, 

(4.3) /ifc+i = M{fik). 

Moreover, {/^fc}fceN+ converges weakly to fi and thus 

(4.4) hm ^k{A) = fi{A) 

k—^oo 

for every set A satisfying fi{dA) = 0. 

Proof, (i) and (ii) follow directly from the construction of the algo- 
rithm. 



(iii) By i^, (|3^, ([3^, and (|3^, it suffices to notice that the 
SSC ensures the existence, for any cylinder set Ei^ of generation 
k, of a unique point x E Ak C] E^^. The existence is clear by 
construction and the uniqueness holds because the cardinality 
of Ak is equal to the number of cylinder sets Ei^. 
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(iv) Let X G Ak+i. By the SSC, there exists a unique j E M such 
that f^^{x) G Ak. Moreover, if ifc+i(x) = {ii{x), ...,ik+i{x)) G 
M^~^^ is the unique sequence of length k + 1 such that x = 



fik+i( x)iy) for some y G Ai, then ii{x) = j and thus, by (2.4) 
and (|310|, 



M(/ifc)(x) = r>fco/-i( 






j^i2{x),...,ifc+i(x) - '^ife+iW - fJ'k+l[X) 



This proves (4.3) 



The weak convergence holds because the sequence {M'^(z/)} 
converges weakly to the invariant measure /i for every compactly 
supported probability measure u (see [3]). In particular, {/ifc} = 
{M'^(/ii)} converges weakly to /x. Finally, this is equivalent to 
(4.4[) (see, for example, [5] Theorem 2.5.11). D 



4.5. Lemma. // {xq, d) E E x [c, 1] is such that 

{2dy 



(4.6) 
thendB{xo,d)nE^^. 



C'E 



fi{B{xo,d))' 



Proof. Suppose that this is not the case. As both dB{xo, d) and E are 
compact sets, there exists < e < dist{E, dB{xo, d)) such that 

(5(xo, d) \ B{xo, d-t))nE = %. 

Thus, iJ,{B{xo,d)) = fi{B{xo,d — e)), contradicting the minimality of 
{xo,d). D 



Observe that for any pair (a;o, d) E Ex [c, 1] satisfying (4.6 ), Lemma 4.5 
guarantees the existence of a point y E E such that 



(4.7) 



C'E 



{2dy 



{2dist{xo,y)y 



l^{B{xo,d)) iJ,{B{xo,dist{xo,y))' 



4.2. Convergence. In this section we show the convergence of the 
algorithm described in section |3} We do it in two steps. First, given 
xo, r, and y as in (4.7), we prove, in Theorem 4.9, the existence of a 



sequence {{xk,dk)}^=i in AkX Dk such that 



ruk := hk{xk,dk) : = 



(24 



-)■ 



{2dy 



IJkiB{xk,dk)) iJ,{B{xo,d)) 



Kixo,d)=C'iE). 
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However, the algorithm's sequence {(xfc, c/^)}^^ is constructed by choos- 
ing on the A;th step a pair {xk, dk) G Ak x Dk such that 

(4.8) 

(2d Y 
rfik := hkixk, dk) := — '' ~ = min{/ifc(xfc, 4) : (x^, 4) e AkXDk}. 
fik{B{xk,dk)) 

Actually, f or e ach k G N^ there might be more than one pair {xk, dk) 



satisfying (4.8). Thus, if we chose a sequence {(a;fc,4)}fcLi from the 



set of pairs selected by the algorithm, it does not need to coincide 



with the sequence {(xfc,4)}fcLi given in Theorem 4.9, Therefore, one 
needs to show that the sequence of minimal values rhk converges to the 
minimum C^lE). In computer experiments we have observed that in 
many cases the sequence hk{xk,dk) is not monotone. Nonetheless, in 



Theorem 4.13 we prove the desired convergence. 



From now on, we use the notation given in (4.8) and write rmax '■ = 
max{ri, ...,rm}. We refer to the sequence rhk as the sequence chosen 
by the algorithm. 



4.9. Theorem. Let xq, r, and y be as in (4--^)- Then, there exist 
sequences {xk} and {yk} such that Xk, yk G Ak for all k G N+, Xk — )■ xq, 
yk -^ y, dk := dist{xk,yk) -> (i = dist{xo,y), and 



lik{B{xk,dk)) fi{B{xo,d))' 
Proof. The existence of the convergent sequences {xk}, {Vk}, and {4} 



is given by Lemma 4.1| (i 



As dk ^ d and, by condition (3.7), /ifc(-B(xfc, 4)) is bounded away 



from zero for all k G N''", we need only to show that 

(4.11) \fik{B{xk,dk)) - fi{B{xo,d))\ ^ as k ^ oo. 

Let 5 > and let (ps be a compactly supported continuous function 
on M" such that < (ps < 1 everywhere, 05 = 1 on B{xo, d + 6), and 
05 = off of B{xo, d+26). By the weak convergence of /i^ to fi (Lemma 



4.1 (iv)), we have that 



bsdfi = hm / (psdfik. 

k-^oo 



Moreover, as ^{dB{xo,d)) = 0, 

\fi{B{xo,d)) - I (j)sdfi\ 
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is small when S is small enough. So (4.11) holds if 

\lik{B{xk,dk)) - j (psdjikl 

goes to zero as 5 — ;■ and A; — > oo. For each 5 > Q take k^ = kQ{5) G N^ 
such that for all k > k^, B{xo, d — 6) C B{xk, dk) C B{xo, d + 6), and 



max 



r'l,^^ < 5. Hence 



(4.12) 

and so 

\fik{B{xk,dk)) 



(psdfik > f^k{B{xk,dk)), 



(j)sd^k 



(psdfik - iJ.k{B{xk, dk)) < 

< fik{B{xo, d + 26)\ B{xk, dk)) < fik{R{26)), 

where R{6) := B{xo, d + 6)\ B{xq, d — 5). 

There remains only to show the convergence of /Zfc(-R(2(5)). Let 
{xi,...,a;j} be the set of points in Ak fl R{25) and, for j = 1,..., J, 
denote by hij) G M'^ the unique sequence of length k satisfying 

fikij)iy) = ^j for some y e Ai, 
and by -Eij.(j) the unique cylinder set of generation k such that Xj G 
Ei^(j). Then ( [3lo| ) and (|42|) imply 

J J 

/^fc(i?(2<5)) 



J 

j=i 



Hk[Xj 



< 



J=l 

< yu({^i, C ^ such that E;^ n i?25 7^ and i^ G M''}) < yu(i?35) 

The last inequality holds because |£^i^| < r'^^.x — ^ fo^ ^11 ifc G -/W^*^ with 
k > ko. This completes the proof of (4.11) as fi{dB{xo,d)) = (see 

m)- □ 

4.13. Theorem. The sequence {'rhk}k&n chosen by the algorithm, given 
in l^.iSj j, converges to C^{E). Moreover, for any /c G N such that 
Hk{B{xk,dk)) = n{B{xk,dk)), there holds 

C\E) < rhk. 

Proof. The second statement in the theorem is immediate, as fik{B{xk, dk)) = 
fi{B{xk,dk)) and (1.6) together imply 

C%E) < hs{xk,dk)) = hk{xk,dk)) = rhk. 

Thus, we need only to prove the convergence. 

For {mfc}fcgN+ as in (4.8), let {rhk^}i^fq+ be a convergent subsequence 
of {Thk}keN+ and m := limj_^oo ""^fci • Since the algorithm only takes 
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values in the set A x (D n [c, 1]) , the sequence {mk}k<^n+ is contained 
in a compact set and it is enough to show that m = C^{E). 

Let A = m — C^^E) and for {mfe-}jgN+, let {mfc-}jgN+ be the cor- 
responding subsequence of the convergent sequence {'mk}ki^^+. Then, 
given < e < |A|/2, there exists io G N"*" such that for any i > io, 



(4.14) 
(4.15) 



\hXxk,,dkJ-C'iE)\ < e/2, 
\hkXS:k,,dkJ -rh\ < e/2. 



First we want to show that m < C'^{E). If this is not the case, then 
A > e > 0. Together with (4.14) and (4.15) this implies 

hk,{xk,,dki) < C\E)^e/2 = m-A+e/2 < hkXxk,,dk,)-X+e < hkXxk,,dk,), 

contradicting the minimality of hkXS^k.jdki)- Therefore, m < C^{E). 

With the aim of showing the reverse inequality, suppose that there 
exists ii G N"^ such that for any i > ii, 



(4.16) 



hkA^k,,dkJ > hsixk,,dkj -e/2. 



and assume that m > C^{E). Then A < — e < and thus, (4.14), 
( |4.15[ ), and ( |4.16[ ) give 

C^(E) = m-A > hk,{xk,,dk,)-e/2-\ > hs{xk,,dk,)-t-\ > hs{xk,,dkj, 

contradict ing t he minimality of hs{xo, d) = C^{ E). T here remains only 
to prove (4.16). Notice that, as dk^ G [c, 1], (2.3) implies that it is 
enough to show that for any sequence {(xfc, dk)}^=i in A x (D fl [c, 1]) 
and (5 > 0, there exists k^ G N"*" such that for any k > ko, 

(4.17) \iik{B{xk,dk)) - fiiB{xk,dk))\ < 6. 

For every k G N+ satisfying that r'^^^ < c, let 



Gk — 
Pk = 

Rk = 

Then, 



{Ei, : ifc G M'^ and E;, C B{xk,dk)}, 

{Ei, : ifc G M^ Ei, n B{xk, dk) ^ and E,^nE\ B{xk, dk) ^ 0}, 

{Ei, ePk-. Akf\ Ei^ n B{xk, dk) ^ }. 



Rk 

Er]B{xk,dk) 

ftkiB{xk,dk)) 

fx{B{xk,dk)) 



C Pk, 



GkU{B{xk,dk)nPk), 
fJ'k{Gk) + fik{Rk), and 
^{Gk)+^{B{xk,dk)nPk). 
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Moreover, by (4.2), jj,{Gk) = Hk{Gk) and jji{Rk) = jJ'kiRk)- AH this 
together with the triangle inequahty, gives 

\Hk{B{xk,dk)) - n{B{xk,dk))\ = \l^k{Gk) + fJ^kiRk) - K^k) - fx{B{xk,dk) n Pk)\ = 

= \fik{Rk) - fj'{B{xk,dk)nPk)\ < 

< fi{Rk) - /i(5(xfc, 4) n Rk) + /i((5(xfe, 4) n Pk) \ Rk) < 

< fi{Rk) + fiiiB{xk, 4) n Pk) , \Rk) < KPk) < 

< fiiB{xk, 4 + rLx) \ Bixk, 4 - rLx)), 

where the last inequality holds because Pk C B{xk, dk + r'^^^^) \ B{xk, dk — ''"max)- 
Thus, we want to show that fi{B{xk, dk+r'^^^)\B {x k, 4— ''"max)) < (^ if ^ 
is big enough. This is true, since E x [c, 1] is a compact set and the con- 
tinuity of fi, together with the fact that fj,{dB{x, d)) = (see pTj), im- 
plies that /i(-B(x, (i+r^j^x)\-S(^5^~''^max)) < SioTajYy{x,d) e Ex[c,l], 



if k is big enough. This proves (4.17) and concludes the proof of the 



theorem. D 



4.18. Remark. Note that the second statement in Theorem |4.13| pro- 
vides an upper bound for C^{E). If on the fcth step the algorithm 
selects a ball containing all the fcth generation cylinder sets that it in- 
tersects, fhk provides us with an upper bound for C^{E). In general, 
the existence of such a ball at step k is not guaranteed as it depends on 
the geometry of the set E. However, fhi will always provide an upper 
bound for C^{E) as the optimal ball selected on the first step already 
contains the whole set. Observe that this value is also an upper bound 
for the spherical Hausdorff measure. Even when we do not have such 
a ball, we can obtain natural upper bounds by estimating the value of 
n{B{xk, dk)), as, for every k eN there holds 

C^E) < ^^^^'^1 
fi{B{xk,dk)) 

(see section [5] for further discussion and examples). 

5. Examples 

The first part of this section is devoted to showing that the algorithm 
recovers easily examples from the literature in which the centered Haus- 
dorff measure has been computed. We have checked how it performs 
for the Cantor type sets analyzed in [1], [21], |22], and [23]. Here we will 
illustrate the accuracy of the algorithm by testing it on three of these 
examples. We will also see that the algorithm suggests the correct min- 
imizing balls for hs{x,d). In the cases where C'^{E) remains unknown 
the algorithm still provides candidates for minimizing balls (and hence 
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for optimal coverings, by [9] and [10]). These potentially minimizing 
balls yields conjectural values for C^i^E) that, in the cases considered 
here, are rigorous upper bounds (see Theorem 4.13, Remark 4.18, Ex- 
ample 5.9, and Conjecture 5.14). In Remark 5.7 we point out that, in 



some cases, we can use the balls selected by the algorithm to obtain 
also the corresponding lower bounds, although the proofs we have so 
far are case specific and are not simple. In the last part of this section 
we present other examples not included in the literature together with 
some conjectures on the precise values of C^i^E) for these cases. 



hJx.r) 



(2r)^ 



Recall that, by Theorem 1.4, our goal is to find the minimal value of 



when X G E and c < r < R. For this the algorithm 



selects, at the kth step, some of the pairs {xk,yk) £ ^fc x ^fc satisfying 



rrik = hk{xk,dk) 



{^dkY 



Hk{B{xk,dk)) 



m\n{hk{xk,dk) : (xfc,4) e AfeXD^}, 



where 4 = \xk - VkV ^o = {^i e M" : fi{xi) = Xi, i = 1, ...,m} is the 
set of fixed points for the similitudes in the IFS, and A^ = S'^{Ak_i) 
(see section 3J). We shall refer to B{xk,dk) as the minimizing ball (or 
interval) . 



1. A-Cantor sets in the real line 

Let K{X) be the attractor of the iterated function system 
{fo{x) = Xx, fi{x) = 1 - A + Aa;}, x e [0, 1]. In |21j it is proved 
that if < A < 1, then 



(5.1) 



C%K{\)) = 2%l-\y, 



where s 



■log A 



is the Hausdorff dimension oi K{X). 



Let A = |. Then K{^) is the middle-third Cantor set. 
The following table gives the results of the algorithm for K{^) 
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after 14 iterations. 





rrik 


{xk,yk) 


dk 


B{xk,dk) 


k = 


1.54856 


(0,1) 
(1,0) 


1 


[-1,1] 


k = l 


1.03238 


V3, 3^ 
fl 2N 


0.333333 


[1,1] 

[0,1] 


2<k<l3 


1.19902 


(1,0) 

(1,1) 


0.666667 


[o,|] 

[-1,1] 



Note that, from (5.1) with A 



|, it follows that 



^^(^(3)) 



, log 2 
4 log 3 



1.199023, 



so the algorithm has found the exact value of the Hausdorff 
centered measure already at the third iteration! A simple com- 
putation shows that, for any k > 2, 



hi 



1 2 
3'3' 



hk 



,2 2. 
3'3' 



(2i)^ 



4« 
~2 



C'{K{\)). 



Notice also that hk{\, |) does not depend on k. In fact, one 
checks easily that 

12 2 2 1 



(5.2 



.3,3/ .V3,3, 

Thus, -B(|, |) = [— |, 1] and -B(|, |) = [0, |] are the actual 
minimizing intervals for hs{x,r) and the algorithm has found 
them already on the second step. 
2. (Ai, A2)-symmetry Cantor sets in the real line 

Let K{Xi,X2) be the symmetric Cantor set defined as the 
attractor of the iterated function system {/o(x) = Aix, /i(x) = 
X2X + ^, /2(x) = 1 - Ai + Aix}, X e [0, 1]. By [4J we know 
that if 

1 - 2Ai - A2 



>A, 



then 



(5.3) 



C^(K(Ai,A2)) = l, 
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where A = max{Xi, A2}. 

Here is the output of the algorithm when \i = I and A2 



1. 

5' 





nik 


{xk^Vk) 


dk 


B{xk,dk) 


k = 


1 


ihD 


0.5 


[0,1] 


k = 1 


1 


a,i) 


0.5 


[0,1] 


k = 2 


1 


(i>i) 


0.5 


[0,1] 


k = 3 


1 


(i>i) 


0.5 


[0,1] 



There is no need to iterate more as this stabihty is empirical 
evidence that the algorithm has found the optimal ball. Indeed, 
we have tested the algorithm with this kind of symmetric Cantor 
sets for different values of Ai and A2 and, in all our trials, [0, 1] 
appears as the minimizing interval for every k > 0. Thus, [0, 1] 
is the natural candidate to be the minimizing ball for hg as well. 
In fact, it is easy to use the exact value of C"*(/C(Ai, A2)) given 



in [1] (see (5.3)) to check that this is the case. To see this, 
notice that for any Ai, A2 as in (5.2), Aq = {0, 1, 1} and 

'2'2^ 



hi 



.1 1, 



h. 



l = C%Ki\^,\2)). 



•2'2' """^2'2' 1 

3. Cantor type sets in the plane 

Let S be the attractor of the iterated function system {/i, /2, /s, f^} 



where fi{x) 
< A, < 



1 



: XiX + bi, i = 1,2,3,4, X = (a;i,X2) G R, 
h = (0,0), h = (1-A2,0), 63 = (I-A3, 1-A3), 



and 64 = (0, 1 — A4). In [23j it is proved that if the parameters 
■^1, -^2, A3, A4, satisfy the conditions 

1 



1 



2Ap< 



n 



(1-Afe)» 

{1,2,3,4} 



< 



) 1-= , where k ^ p, k, p G {1, 2, 3, 4}, 
where k ^ i, k ^ j , i j^ j, k, i, j E 



then 



(5.4) 



C'iS) = D 



-1 
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where 
(5.5) Dmax = max 






i<*<4i (2^2(1 -Ai))^ 
v^Ai 



mm 



^k + K + ^ 



(min{l-A„l-A,})^'(v^(l-Ap))« 
{k,i,j,p) G {(1,2,3,4), (2,1,3,4), (3,2,4,1), (4,1,3,2)} 

and 

< dimH(^) = s < 1. 
Next, we test the algorithm for the case Ai = A3 = ^ and 
A2 = A4 = ^. In this case we have that s = log2o("\/3 + 1) ~ 



0.335495 and 



C'{S)= y^] ~ 1.3932 



see [23] . 

The next table shows the results given by the algorithm after 
six iterations. 

Table 1. 
Output of the algorithm for the case Ai = A3 = ^ and 







A2 — A4 — ^. 






k 


ruk 


{xk^Vk) 


dk 


B{xk,dk) 


k = 


1.4174 


((0,0), (1,1)) 
((1,0), (0,1)) 


v^ = 1.41421 


B{{0,0),V2) 
((1,0), (0,1)) 


l<k<5 


1.39321 


((|j,i),(l>0)) 
((i,|j),(0,l)) 


19^ = 1.3435 


B(( 1 19) 19a/2x 
VV20' 20^' 20 ' 

^//19 IN 19v^N 
VV2O' 20^*' 20 '' 



We note that the pattern of the previous examples is re- 
peated here. There exists ko such that, for all k > ko, the 
same balls are selected. The results in Table [T] show that for all 
k > 1, the optimal balls for the discrete density function h^ are 
B{{^, I), 1|^) and 5((1§, ^), ^f), so these should be con- 
sidered as the natural candidates to minimize hg. As before, by 
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fact the minimizing balls: 



hs{{ 



(|5.4|), a simple computation shows that these candidates are in 

19 1 , 19^2 _ 

19 1 , 19v^, 



1 19 19^2 
20' 20'*' 20 ' 



(5.6) 



hsii 



20 20' 
1 19. 



20 ' 
19^2, 



20 20' 
I9V2' 

"To" 



20 ' 

C'{S). 



hkii 



20' 20^' 20 



The third equality in (5.6 ) holds because fik{B{{^ — "> ^^^^ 

A''fc(-°((20' 20''' "20"'''' 



^20' 20^' 20 



)) 



1 for any A; > 1. 



h 



C'E 



r i2ry 
min < —, — ; r- 

(2(1 -A)r 



X E E and c < r < R 



M5(l-A,l-A)) 
The proof reduces to showing that 

{2ry 



2'(i-xy 



(5. 



/i(5(A,r)) 



> 2^(1 - A)' for any r G [1 - 2A, 1 - A], 



' 20 



5.7. Remark. As we have seen in the preceeding examples the algorithm 
is also helpful for finding natural candidates for optimal balls. For 
instance, an output repeated for all iterations k > k^ is a. natural 
candidate. Once the algorithm has suggested such an optimal ball, 
this can be used to provide new proofs of (5.1), (5.3), and (5.5). We 



explain the main idea for the first case. The other cases can be treated 
in a similar fashion. 

A-Cantor sets in the real line. We noticed that when we applied 
the algorithm to different K{X) we obtained two values for {xk,yk), 
namely (A, 1 — A) and (1 — A, 1 — A). Once we have a candidate for 



the optimal ball, we can recover (5.1 ) by showing that the minimum of 
X, r) is attained for the pairs (A, 1 — A) and (1 — A, 1 — A), i.e. 

(2(1 -A)r 



/i(5(A,l-A)) 



since the other cases are trivial. The inequality (5.8) can be proved by 
hand using the geometry of the set. 

For the (Ai, A2)-symmetry Cantor sets in the real line the same 
method can be used to show that the optimal interval is actually the 
one chosen by the algorithm, namely the interval [0, 1]. Finally, for the 
attractor S of the third example, the candidate for an optimal ball is 
5((Ai, Ai), ^2(1 - Ai)) whenever Ai > A2 > A3 > A4. 
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We note that the results in [1], [21], [22], and [23] are obtained us- 
ing the relation between the upper spherical density of the invariant 
measure and the centered Hausdorff measure. More precisely, in these 
references it is shown, that from the fact that, if C^{E) < oo, there 
holds 

d'{C^E,x) = hmsup ' ' = 1, 

for almost all x & E (see [IS]), there follows 

C^(E) = (J^(/i,x))-\ 

for almost all x G E, where /i = ^J^. is the invariant measure. So 
they proved that 

cf (ii, x) = lim sup — — = — ; -^ , 

for every point x in a certain set of full /z measure. Notice that the 
method suggested by the algorithm is different since it is not necessary 
to pass to the limit. 

Next, we present the results obtained by the algorithm for a family 
of Sierpinski gaskets whose centered Hausdorff measure is not known. 

5.9. Example (Sierpinski gasket S{r)). Let S{r) be the attractor of 
the system ^ = {/q, /i, /s}, where 

(5.10) foix,y) = r{x,y), 

(5.11) fiix,y) = r(x,|/) + (l-r,0), and 

(5.12) f2{x,y) = r(x,|/) + (^(l-r),(l-r)^)). 

After testing many examples, we have discovered that for r < 0.25 the 
algorithm quickly stabilizes at [2(1 — r)(r^ +r + 1)2]*, where s = ~^°^ 
is the Hausdorff dimension of S{r). We conjecture that, in these cases, 

C'{S{r)) = [2(1 - r){r^ + r + 1)^]^ 



We illustrate this with Table |5.9[ which shows the output of the algo- 
rithm for r = 0.2. 



algorithm for centered hausdorff measure 
Table 2. 

Output of the algorithm for the case r = 0.2. 
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k 


mfe 


{xk,yk) 


dk 


B{xk,dk) 


k = 


1.60504 


((^,f),(l,0)) 
((^,f),(o,o)) 


1 




k= 1 


1.51231 


((Ii'0),(^,f)) 


0.916515 


5((|)>0),^) 


2< A;< 8 


1.48326 


((1. yl) (A. yl)) 

VV5O' 50 >'' VlO' 2 ^^ 


0.890842 


5((|),fX^) 



The pattern observed in our trials for these Sierpinski gaskets is the 
following. For any k > 2 the optimal ball chosen by the algorithm 
is B(x, d{x, X2)) = B(x, (1 — r)\/r'^ + r + 1), where x = /o(/i(^2)) and 
X2 = f{x2) = f{x2-, 1/2) € ^ is the fixed point of /2(a;, y) (see FigurejTj). 
Observe that S{r) C B{x,d{x,X2))- Hence, for any fc G N, 

jJk{B{x, d{x, X2))) = /i(-B(S, d{x, X2))) = 1. 



Then, by Theorem 4.18, we have a rigorous upper bound for C'^{S{r)), 
namely 



C'iSir)) < [2{l-r){r^ + r + l)-^Y 



Vr <0.25 



5.13. Remark. Note that in all previous examples the minimizing ball 
(or interval) is big enough to cover the whole set. That is, if B is the 
optimal ball, then fi{B) = 1. The reason for this lies in the fact that 
the contraction factors are very small compared with the distances that 
separate the cylinder sets. 



In the case of the Sierpinski gaskets S{r) (Example 5.9), if r < 0.25 



then the candidate for an optimal ball is selected already at the second 



step because r is still quite small in comparison with c (see Theorem 1.4 



for notation). With larger values for r, the centers of the selected balls 
are the same after the second step, but the end point varies slightly. 
Therefore, the value of m^ changes slightly after the second step, but 
in a few steps the output is stable to three decimal places. This leads, 
for example, to the following conjectures for the Sierpinski gasket S'(|) 
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s= -logS/logr 




Figure 1: 
Minimizing ball selected by the algorithm for S{r) with r < 0.25. 



and the well-known planar ^-Cantor set, C(|), both self-similar sets 
with dimension 1. 
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5.14. Conjecture. Let S'(^) be the attractor of the system \E' = {/o,/i,/2} 



^3 

given in (5.10) with r = ^ and let C(^) be the attractor of the system 



* = {90,91,92,93} where 



9o{x,y) = 


= ^ix,y), 


9iix,y) = 


- \i^,y) + il,o) 


92{x,y) = 


= i(x,y) + (0,^) 


93{x,y) = 


= i(^'^) + 4^; 



The conjectural values given by the algorithm for these self-similar sets 
are 

(1) C\S{l/3)) ~ 1.537. 

(2) Ci(C(i))^1.95. 

We will analyze the planar ^-Cantor set to show how to obtain a rig- 
orous upper bound from the results given by the algorithm. The Sier- 
pinski gasket >S'(|) can be treated in a similar manner. 

First we present a table showing the results of seven iterations of 
the algorithm (see Table |3]). For simplicity we focus on only one of the 
minimizing balls obtained. In this case we believe that the description 
of the pairs {xk,yk) by their corresponding codes is more illustrative 
than the description by their coordinates, so we will use the following 
standard notation: for every k eN and iii2.--ik G M^, 



Xiii2...ik — 9iii2—ik-i\-^ik) ~ 9ii ° 9i2 ° ••• ° 9ik-l\-^ik), 



where Xi^ = 9iki^ik), ^ = {0,1,2,3}, and M'' = {iii2...ik ■ ij e 
M ^ j = l^...,k}. Observe that if it = it+i = ■■■ = ik for some 
< t < /c, then a;j^J2...jj, = Xj^J2...j^. We shall always use the shorter 
notation. 
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Table 3. 



Output of the algorithm for C(^). 



k 


rrik 


{xk,yk) 


dk 


k = 


2.66667 


(a^ca^s) 


1 


k = 1 


1.92296 


(3^20, a^os) 


0.901388 


k = 2 


1.95814 


(x2o, 3^020) 


0.795495 


k = 3 


1.95542 


(a;2o,a;3) 


0.790569 


k = A 


1.95306 


(a;2o,a;3) 


0.790569 


k = 5 


1.95388 


(X20, 3^020023) 


0.790662 


k = 6 


1.95417 


(X20, 3^020021) 


0.790662 



For simphcity, we are going to use the value of ma in Theorem 4.13 to 



4.13 



we have 



obtain the upper bound for C^(C(|)). It easily seen that fi^^B^xs, d^)) 

fj,{B{x2,,d3)) = I + p- + |j (see Figure 2), so by Theorem 
that 

C^'(C^(^))</3(^3,^3) = 1.95542. 
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/ "1 rl h r T n n r 

1 3 d b c 








Figure 2: 
B{x4, (^4): Minimizing ball selected by the algorithm at step k = A for 
C(|). 
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